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Abstract 

Given the spectrum of a Hamiltonian, a methodology is developed which employs the Landau- 
Ginsburg method for characterizing phase transitions in infinite systems to identify phase transition 
remnants in finite fermion systems. As a first application of our appproach we discuss pairing in 
finite nuclei. 
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Recently it has been pointed out that empirical evidence exists for a pairing phase tran- 
sition to occur in symmetric nuclear matter at normal nuclear densities at T c ^8 Mev 
Here the energy density and the specific heat have been obtained from a finte tempera- 
ture extension of the semi-empirical mass formula e 3]. A Landau Ginzburg treatment of this 
transition together with a simple pairing calculation strongly suggests that its origin is due 
to the existence of a paired superconducting phase at temperatures below .8 Mevj^J. This 
result is not surprising as finite temperature BCS calculations in nuclear and neutron matter 

have suggested that such a phase should exist 0, ?l A Is > H HQ, HQ HQ- 
one of the intriguing questions that arises is whether or not a remnant of this phase 
transition survives in finite nuclei. Clearly a universal feature of finite nuclei is a significant 
change in the density of states at excitation energies of 10 MeV or lessjl?|. At lower exci- 
tation energies the spectrum of most nuclei is sparse and dominated by a relatively small 
number of collective states. With increasing excitation energy, the independent particle 
degrees of freedom dominate and the density of states grows exponentially. As the mass 
number increases, the low-lying collective portion of the energy spectrum becomes more 
compressed and an abrupt change in the many particle density of states occurs at lower ex- 
citation energies. It has, therefore, been suggested that a collective to non-collective phase 



transition occurs in finite nuclei 



181 ] . Strictly speaking it is incorrect to speak of phase 



transitions in finite systems. No one can deny, however, that transitions between different 
regimes do take place in these systems, and that they are more or less abrupt. (A similar 
situation is known in biophysics. The helix-coil transition in certain biological molecules in 
solution occurs with the temperature width of around 5 degrees.) This issue, particularly 
in deformed systems, has been clouded by the fact that finite temperature mean field calcu- 
ations have suggested that this phase transition is simply due to a drastic change of shape 



[19]. 



The deformed-to-spherical 
exact canonical calculations 



shape transition seen in these calculations is not seen in 
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22iQ and may be an artifact of the finite temperature 



mean field approximation and also depends on the volume of the system [24, 25]. In spite of 
the fact that the canonical partition function above the critical temperature is dominated 
by the single particle degrees of freedom a few collective states still contribute and are ex- 
tremely important in calculation of shape dependent parameters. Recent calculations of the 
ensemble average of the quadrupole moment squared Q' 2 '- indicate that it is discontin- 
uous in the finite temperature mean field approximation, while no discontinuity is observed 
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in the canonical calculations j^. In both cases this quantity does not appear to vanish at 
the critical temperature. It should also be noted, however, that when thermal fluctuations 
in the shape dependent order parameters are taken into account, either by macroscopic or 
microscopic procedures, reasonable agreement with the exact canonical calculations 2(| is 
obtained. With increasing temperature, however, it is expected that these collective degrees 
of freedom will eventually completely dissolve, presumably below the critical temperature 
of the liquid-to-gas phase transition. 

However, it should be noted that model studies in a SU(2) xSU(2) system show that, in 



the thermodynamic limit, 
of a true phase transition 



;his system exhibits a singularity in the specific heat characteristic 
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281 ] . Furthermore, the remnant of this singularity remains in 
the form of a peak in finite sytems of this type. The presence of this peak has been used to 
map out the phase structure in such a model 



In spite of the fact that in many microscopic variational calculations (see for example 30]) 
pairing transitions appear to take place they are difficult to identify in exact shell model 
calculations [3Jj. We propose in the present work to identify the existence of a pairing 
transition empirically in finite nuclei in the following manner. In the Canonical Ensemble 
the partition function for a nucleus of mass A is given by 

Z(A,T) = J> exp{-(3E i ) + / dE g AZ (E) exp{-(3E) (1) 

where (3 — ™, <?j = 2jj + 1 is the spin degeneracy factor, Et the energy of the ith state of 
the nucleus and gA,z{E) its level density with 

1^ at 

where U = E- P(Z) - P(N) and a = A[0.00917S + 0.142] per MeV for undeformed nuclei, 
or a = A[0.00917S' + 0.120] per MeV for deformed nuclei. A is the mass number and S is 
the spin energy. Given experiment information about the bound states of a nucleus and an 
experimental fit to its continuum level density, the partition function of a nucleus can be 
determined empirically 

From Z it is easy to determine as a function of T the excitation energy 

E = -J^nZ (3) 

and the specific heat 



Cv = §. (4) 

A second order phase transition such as a pairing phase transition in an infinite system 
should lead to a discontinuity in Cy at T c . Clearly such behavior is not posible in finite 
sytems. Only a remnant remains in the form of peaked structure in the specific heat. In 
order to ascertain whether this structure is consistent with a pairing phase transition we 
make use of Landau- Ginzburg theory to demonstrate that a description of the condensed 
or paired phase close to the critical temperature, T c , from information in the normal or 
uncondensed phase obtained from experimental data can be can be accomplished. The 
magnitude of the remnant of this discontinuity in Cy is compared with a simple analytical 
calculation for a pairing phase transition jj] corrected for finite size effects j^]. 

Landau and Ginzburg have provided a simple theory of phase transitions which approx- 
imates the free energy in the region around T c and is most useful in analyzing the thermo- 
dynamics in this region. In particular, using only knowledge about the uncondensed phase 
one is able to make predictions about quantities in the condensed phase,such as specific 
heat, magnetic susceptibility and compressability. Moreover, Landau- Ginzburg theory can 
be derived from microscopic considerations ^J. 

In the Landau Ginzburg formulation it is necessary first to determine an expression for 
the free energy F(T) in both phases. In the following the subscript 1 will refer to the lower 
temperature (condensed) phase, and 2 to the higher temperature (uncondensed or normal) 
phase. In the uncondensed phase, a quadratic form for the energy, which follows from a low 
temperature Fermi gas, is used approximation of a normal Fermi liquid, 

E 2 (T) = a 2 + k 2 T\ (5) 

where a 2 and k 2 are constants. From the relations for the specific heat terms of E and the 
entropy S, 

dE _ OS 

Q y-of- T of ' (6) 

one obtains the entropy in the uncondensed phase, 

S 2 (T) = C 2 + 2k 2 T, (7) 
where C 2 is an unknown integration constant which later cancels out of the calculation. 
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From eqs. (J^J) and (J7J) the free energy in the uncondensed phase is given by 

F 2 (T) = a 2 -C 2 T-k 2 T 2 . (8) 

pi 

The free energy in the condensed phase is obtained from the Landau expansion 4] for the 
free energy in terms of an order parameter rj which goes to zero at the transition to the 
uncondensed phase. This order parameter vanishes at a critical temperature T c . The free 
energy expansion to order r\ is 

F 1 (T,? ? ) = F 2 + Arf + Brf. (9) 

Here A and B are functions of temperature and it is assumed that the states with rj=0 and 
rj t^O are of different symmetry. In this case it can be shown the linear term in r\ must be 
set equal to zero and if the critical point is also a stable point, e.g. if Fi as a function of r/ is 
a minimum at rj =0, then the third order term in 7] should be zero and at the critical point 

D 

A = B > 

The order parameter is determined by requiring the condensed phase to be stable below T c 
(i.e. Fi should be minimized w.r.t. rj). This leads to 

F, = F 2 - £ (10) 

Furthermore, since A is of opposite sign in the condensed and uncondensed phases, while B 
is strictly positive (^J , the lowest order expansion of A in T — T c can be parametrized as 



A(T) = a(T-T c )2^/B(T c ). (11) 

Note especially that a > is an essential requirement following from the phase diagram Q]. 
Substituting for A(T), the free energy per nucleon near T c is given by 

Fi(T) = (a 2 - a 2 T 2 c ) + (2a 2 T c - C 2 )T - (k 2 + a 2 )T 2 , (12) 

where F 2 is taken from eq. (JBJ). 

From eq. ()12j) . the energy in the condensed is easily determined to be phase near T c , 

Ei(T) = (a 2 - a 2 T 2 ) + (a 2 + k 2 )T 2 (13) 

= a x + k{T 2 . (14) 



Comparing this to the uncondensed phase (eq. (j3J)) we note that the T dependence is also 
quadratic, but has a larger coefficient. Thus the specific heat is discontinuous across the 
phase transition, and is necessarily larger (ki > k?) in the condensed phase. 

We now apply our theoretical results to the following even-even nuclei, namely, 20 Ne, 48 Ca, 
88 Sr and 208 Pb. These nuclei span a large portion of the mass spectrum and have a reasonably 
well determined energy spectrum both in terms of their energies and the corresponding 
angular momentum assignments. For each of these nuclei, the partition function (see eq. 
(1)) contains a discrete as well as a continuum contribution. 

The discrete portion of the partition function is determined empirically from the mea- 



sured energy levels in each nuclei jl7|]. At lower excitation energies the angular momemtum 
assignment to each state in the spectrum is generally unique. As the energy rises, pre- 
dominantly near the onset of the continuum, states occur which have an uncertainty in the 
angular momentum assignment. In such cases, we have taken the lowest suggested value of 
the angular momentum. The continuum contribution to the partition function is given by 
equation (J2J) and the value of the parameters used for the even-even-nuclei considered are 
given in Table 1. 

The energy at which the continuum is to be attached to the discrete por tion of the energy 
spectrum is determined via the prescription of Gilbert and Cameron [l7j. This matching 
point was ascertained graphically by first plotting the energy as a function of the number 
of levels. There are two parts in each such graph, viz, the curve corresponding to the low- 
energy bound states, and the curve corresponding to the high energy continuum states. The 
point of tangency (i.e. where the slope is the same in both curves) was then determined 
and selected as the matching point. They found that this could be parameterized by the 
following 

E C = U X + E P (15) 

In the above equation, E c is the energy at the matching point and E p is the pairing energy 
associated with the nucleus under consideration. U x is an additional energy term which 
is found graphically from the tangency point (minus the pairing energy) as a function of 
mass number (see Gilbert and Cameron 3]). From the graphical results one can see that 
there is an upper and lower limit to U x , that the curves are hyperbolas, and that the overall 
behaviour is a decrease in U x with increasing A. The range for U x is listed below. On 
average, 



TABLE I: Parameters used for nuclei 



Nuclei 


P(Z) (MeV) 


P(N) (MeV) 


S(Z) (MeV) 


S(N) (MeV) 


a (MeV" 1 ) 


E c (MeV) 


208pk 


0.83 


0.38 


-8.86 


-3.16 


6.61 


6.25 


88 Sr 


1.24 


0.93 


-16.41 


12.88 


9.65 


7.14 


48 Ca 


1.83 


1.3 


-12.07 


12.13 


6.80 


10 


20 Ne 


2.5 


2.5 


-0.811 


5.633 


4.58 


17.7 



U x = 2.5 + ±f , 

the upper limit is , U x = 2.7 + 

and the lower limit is, £4 = 2.1 + ^p. 
We have used the upper limit in all our calculations, except for 208 Pb as it produced the 
best fit to a straight line in the uncondensed portion of the specific heat (see figures 1-4). 

For 208 Pb it is extremely difficult to determine the density of states. It has been pointed 
out by Gilbert and Cameron that not all nuclei near closed shells could be fitted (i.e. finding 
the tangency point) in the aforementioned manner. For those that could be fitted, a value 
for U x was found that was much larger than the one predicted. 208 Pb is such a case. One can 
determine U x from the graphical results given in Gilbert and Cameron (this value is roughly 
4.8 MeV). Adding the pairing energy to this yields a matching point value of roughly 6.5 
MeV. We used a slightly lower value. The parameters used in the determination of the 
partition function as a function of temperature for the nuclei under consideration are given 
in Table 1. 

For all four nuclei (see Figures 1-4), a well defined peak is observed in the specific heat ( 
equation^) as a function of temperature and the Landau- Ginzburg fits to the specific heat 
are consistent with that of a remnant of a second order phase transition. The discontinuity 
in the specific heat expected in an infinite system appears to be smoothed over at the critical 
temperature due to finite size effects. Furthermore, not unexpectedly at higher tempertures 
above the critical temperature, the specifc heat is linear as a function of temperature. 

Varying the values of the matching points for the four nuclei considered gives rise no 
significant change in quantities like the critical temperature (variation of 0.01 MeV), or to 
the slope of the straight line fit (in the uncondensed phase) to the specific heat at higher 
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Cv 




FIG. 1: Specific heat of 208 Pb as a function of temperature for the parameters given in Table 1. 
The dashed curve is the Landau-Ginzburg fit to the specific heat, whereas the solid curve is the 
specific heat determined from the experimental data. . kl = 10.83 and k2 = 6.75 




FIG. 2: Specific heat of Sr as a function of temperature for the parameters given in Table 1. The 
dashed curve is the Landau-Ginzburg fit to the specific heat, whereas the solid curve is the specific 
heat determined from the experimetnal data. . kl = 22.94 and k2 = 10.79 

temperatures. The numerical value of the observed peak in the specific heat for each nucleus 
does vary slightly. 

In an infinite system for a pairing phase transition, the discontinuity in the specific heat 
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Cv 




T Mev 



FIG. 3: Specific heat of Ca as a function of temperature for the parameters given in Table 1. 
The dashed curve is the Landau-Ginzburg fit to the specific heat, whereas the solid curve is the 
specific heat determined from the experimental data. . kl = 17.10 and k2 = 6.80 




T Mev 



FIG. 4: Specific heat of 20 Ne as a function of temperature for the parameters given in Table 1. 
The dashed curve is the Landau-Ginzburg fit to the specific heat, whereas the solid curve is the 
specific heat determined from the experimental data. . kl = 6.53 and k2 = 3.74 



at T c can easily be analytically determined |4j and is given by 

7C(3)?? 



C S (T C ) - C n (T c ) 



C n (T c ) 



12 
703) 



(16) 
(17) 



= 1.43 



(18) 



and one can define 



A L = C S (T C ) - C n (T c ) 
= 1.43C n (T c ) 



(19) 
(20) 



Clearly in order to calculate the remnant of this discontinuity in C for a nucleus, this result 
should be corrected for finite size effects. A simple way of modeling this is to assume that 
the equilibrated system is contained in a finite sized volume. The sum over the momentum 
states in quantities like the partition function can be approximated then by an integral over 



the density of states in the form 

S„ — 



32 



33] 



p ap± ~pdp H -dp 



2n 2 h 3 



8Trh z 



8nh 



(21) 



where S and L are the surface area and linear dimension of the system of volume V. The 
± correspond to either choosing Dirchlet (-) or von Neumann (+) boundary conditions. In 
the following we shall use Dirchlet boundary conditions and neglect the contribution from 
the linear term. 

From Liftshitz and Pitaevskii ^| it easy to see that there is no change in the numerator 
of equations (16 and 17) when finite size effects are included and that only the surface 
term must be taken into account in the calculation of C n in the denominator. The surface 
contribution to C n is given by 



where 



Csurf 



2S 



a p surf 
_ n 

dT 



LP d P 



i 



fo2 

eXp 2m +1 

oo doc 
z^ 1 exp x +1 



2SmT 2 

8ivh 2 
2SmT 2 

8ivh 2 

2SmT 2 e 2 f r 2tt 2 

[1 + 4 



Fo(z) 



8nh 2 2T 2 ' 



(22) 

(23) 

(24) 
(25) 
(26) 



where e t is the Fermi energy z = is the fugacity and F2(z) is the well known Fermi-Dirac 
integral[32|. From this it follows that the surface contribution to the specific heat for finite 
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nuclei is given by 

C sur f „ VmpjT Snh 

^ ~ 3fr 3 [ AVp f --- 1 {Zf) 

and that 

Cf" te (T c ) - Cfi mte (T c ) V 



r finite /rp\ T/ wp/T c r 1 3SH 



1.43 



[1 



4Vp f i 



(28) 
(29) 



and in same manner one can define 



A{ mite = Cf mte (T c ) -C£ mte (T c ) (30) 
1 4S 

= T-MjCf 6 ^) (31) 
Here the the Fermi momentum, pp, for a nucleus with ^4 nucleons is determined from 

^ " 3(2^)3 (32) 

and for simplicity V = and S = AttR 2 where R = r A^. 

After the Landau-Ginburg fits to the empirical values of the specific heat have been 
obtained (see figures 1-4) it is easy to determine the value of the remnant of the discontinuity 
(hereafter we shall refer to this as the discontinuity) in the specific heat at T c for each of 
the nuclei (see Tablell). As one goes to heavier nuclei T c falls off hyperbolically and reaches 
for 208 Pb a value which is slightly less than that in symmetric nuclear matter (2|. If one 
now assumes that this discontinuity is due to the remnant of a pairing phase transition, 
its value with(A{ mite ) and without finite size corrections (A^) can easily be obtained from 
from equations (|3T|) and fl2()l) . In spite of the simplicity of the pairing calculation(the density 
of states and the pairing interaction are assumed to be constant) and the difficulties with 
the continuum contribution (for 208 Pb) the value of the discontinuity calculated with or 
without finite size corrections is in reasonable agreement with the empirical value for the 
three heavier nuclei. Note that in the heavier nuclei the finite size corrections are not large. 
The discrepancy in all cases is less than w 40%. Only in the case of 20 Ne is the discrepancy 
about 60%, more than 1.5 the value for heavier nuclei. This strongly suggests the existence 
of a pairing transition in the heavier nuclei which is not present in 20 Ne. In the latter case 
this evidence for a second order transition which probably is shape related Q]. 
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TABLE II: The values of the discontinuity in the specific heat atT c 



Nuclei 


A 


A L 


^ finite 


T c (MeV) 


208p b 


7.15 


10.22 


10.27 


0.66 


88 Sr 


22.67 


14.81 


14.92 


0.72 


48 Ca 


21.59 


18.02 


18.27 


1.0 


20 Ne 


12.95 


20.88 


21.46 


2.11 



From the spectra of four even-even nuclei 20 iVe, 48 Ca, 88 Sr and 208 Pb we have empirically 
constructed in the Canonical Ensemble their partiton functions from available experimental 
data and determined their specific heat as a function of temperature . For each nuclei, the 
specific heat displays a prominent peak which may be the remnant of a phase transition. A 
Landau Ginzburg treatment shows unambiguously that this interpretation is not inconsistent 
if the the phase transition is second order. A simple pairing calculation of the magnitude of 
the observed discontinuity is consistent with that obtained empirically for the three heavier 
nuclei. This suggests that a pairing transition takes place in these nuclei. In the case of 
20 Ne such is not the case and the transition may be shape related. 

Lastly we wish to point out that the methodology that we have employed here can be used 
for any finite fermion system. Given the exact spectrum of the system, the Landau-Ginsburg 
method can be utilized to identify the existence of the remnant of a phase transition. In 
cases such as pairing phase transitions, simple analytical calculations can be used to identify 
the nature of the phase transition. 

AP acknowledges support from the Argentine National Science Council. HGM acknowl- 
edges the hospitality of the Physics Department of SUNY at Buffalo where part of this work 
was undertaken. 
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